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Abstract 

Modified Volterra lattice admits two vector generalizations. One of them is studied 
for the first time. The zero curvature representations, Backlund transformations, non- 
linear superposition principle and the simplest explicit solutions of soliton and breather 
type are presented for both vector lattices. The relations with some other integrable 
equations are established. 

Key words: Volterra lattice, Darboux transformation, nonlinear superposition prin- 
ciple, zero curvature representation, symmetry. 

1 Introduction 

Vector equations are an important and rather well studied class of integrable system. 
Among others, we mention but a few works in this area [HOE] containing the examples 
and classification results for the vectorial systems of derivative nonlinear Schrodinger type 
which are in some relation to the theme of our paper. There are also several interesting 
results for the vector differential-difference equations, or lattices, see e.g. [HE], but this 
field seems less investigated. The aim of our work is the study of the vector lattices 

V n , x = 2<K, Vn+i - - (K, V n )(V n+1 - (l) 

V n ,x = {V n ,V n )(V n+1 - F„_i), (2) 

which define two integrable generalizations of the very well known modified Volterra lattice. 
Equation ([T]) was introduced in [B] among the other examples of the multi-component 
lattices related to Jordan algebraic structures. Second lattice is considered here for the 
first time, up to our knowledge, despite of its more simple form. 

The main tool in the study of a nonlinear integrable equation is its representation 
as the compatibility condition for auxiliary linear systems. In the differential-difference 
setting this method was developed in the classical papers [7J [8]. In our paper we restrict 
ourselves by the version of dressing method based on Darboux-Backlund transformations 
and their nonlinear superposition principle. The main results for the scalar lattice are 
given in Section [2j The main body of the paper, Sections [3l [U contains generalizations of 
this method for both vector lattices ([T]), (f2]), as well as the simplest explicit solutions of 
soliton and breather type. 



*L.D. Landau Institute for Theoretical Physics, Chernogolovka, Russia. E-mail: adler@itp.ac.ru 
^Sochi Branch of Peoples' Friendship University of Russia, Sochi, Russia. 
E-mail: postnikovvv@rambler.ru 



I 



A characteristic feature of integrability is the consistency of the equation with an 
infinite hierarchy of other equations. In particular, Backhand transformations define the 
discrete part of this hierarchy and lead to the discrete equations on the square grid. Usually 
one starts this way from the continuous equations of KdV type, however an understanding 
appeared recently that the lattice equations of Volterra type lead to the same result as 
well [TP] lllj . This relation has not been observed in the vector case yet, although the 
discrete equation related to the lattice ([T]) has been introduced in the paper [12], see also 

The continuous part of the picture (Section [5]) is more traditional. It was observed in 
works of Levi [16] and Shabat, Yamilov [T7] that integrable Volterra type lattices define a 
special kind of Backhand transformations for equations of nonlinear Schrodinger type. This 
remains valid for the vector analogs as well. The connection with a two-dimensional lattice 
relative to the Volterra lattice introduced by Mikhailov [18] is of interest, too. Finally, it 
should be noted that the approach based on the continuous symmetries is the most effective 
one in the classification problem of integrable equations, both continuous and discrete [19] . 
The complete classification of scalar Volterra type lattices was obtained by Yamilov [20] by 
use of the symmetry approach, see also |21^ I22j. Some progress in classification of vector 
equations and lattices has been achieved recently [3] [231 12H US]- We discuss some open 
problems in this field in the concluding Section [6] 

2 Scalar case 

2.1 Zero curvature representations 

The following notations for the auxiliary linear equations are used throughout the paper: 
* n+1 = L n * n , V niX = A n ^ n , § n = M n V n . (3) 
Modified Volterra lattice 

Vn,x = (Vn + a){v n+ i - V n -i) (4) 

is equivalent to the compatibility condition L nx = A n+ \L n — L n A n with the matrices 
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The Darboux-Backlund transformation is defined by the matrix 



1 h{a 2 -^\ 2 )-a^X 2 -^)f 2 n -(a 2 -^)Xf n \ 

U « + V 2 ti \ (a 2 - ^)Xf n a^(X 2 - /i 2 ) - /i(a 2 - ^ 2 A 2 )/ n 2 J ' 

Moreover, the compatibility condition L n M n = M n+ iL n is equivalent to the pair of discrete 
Riccati equation for the variable f n : 

1 + fnfn+1 1 + fnfn+1 
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Figure 1. Breather of the lattice ((H; a = 1, 7^ = 7^ = 0.5 + 1.5i, 

fc(D = fc(2) = 1. 



and the condition M„ )X = A n M n — M n A n completes this system with the continuous 
Riccati equation 

fn,x = (-Vn-1 + A^n) fn + (~2 ~ A^) /« + A">n-1 + ^ v n- (8) 

Notice also that the variable f n satisfies, in virtue of equations ©, ([8]), the lattice 

(/i 2 + a/2)(a + / x 2 ^)(/n + i-/„-i) 



fn,a 



(9) 



/i 2 (l + /rH-l/n)(l + /n/»-l) 

Starting from a known solution v n of the lattice (j3|) the common solution of the first 
equation (JTj) and equation (JSJ) is constructed by the formula f n = (pn/^n where \& = (4>, <p) 
is a particular solution of two first equations ^ at A = fi. Then the second equation ([7]) 
defines the new solution v n . 

For example, in order to construct solutions of soliton type one takes v n = 1 as the seed 
solution (obviously, the choice of another constant is equivalent to scaling of parameter a; 
some generalization can be achieved via dressing of blinking solution V2 n = « ; ^2n.+i = /?)• 
The eigenvalues of the matrix L n \ Vn= \ : \ = ^ are defined from the equations 



7i+72 = M + aM 7i72 = l + a 
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and the corresponding solution of the linear equations is 

(we do not consider the case of multiple roots 71 = 72 which leads to rational in n, x 
solutions). The ratio f n = (pn/^n defines the solution of the lattice ^ of kink type 
(provided 71/72 > 0, k > 0) and the substitution into the second equation (|7|) gives the 
soliton of the lattice (HJ). The construction of iV-soliton solution uses the set of particular 

solutions (<j)n ] ifrP) corresponding to the values of parameters j = 1,...,N. 

If a > then the lattice admits the breather solutions corresponding to the pairs of 
complex conjugated points in the discrete spectrum (/Z 1 ) = p^ 2 \ k^ = k^). 



2.2 Nonlinear superposition principle 

The direct recomputing of the variables / is a more convenient way to iterate the Darboux 
transformation than applying the matrices M and recomputing the wave functions. This 
leads to the nonlinear superposition principle of Darboux transformations in the form of 
some Yang-Baxter mapping [26J. Let the variables /A be constructed from the particular 
solutions of the linear systems at \i = and let jf^'iiii— J«) (\ eno i e the variables obtained 
from f$ by consequent application of Darboux transforms with parameters ■ ■ ■ , fj,^ s \ 
Then the permutability of Darboux transformations is equivalent to the following equality 
for the matrices of the form ((6|): 

M(f^\^)M{f^\^ {k) ) = M(fi k ^\^)M(f^\^) 

where a stands for a tail sequence of distinct indices. This equation is uniquely solvable 
with respect to fn' k ' a \ frt'^ and thus the mapping is defined 

M{J, 9: V, v) = 3- — ■ —2 Y\ 2? 7T~- ( iU ) 

V-v (M5 -v})g + av[pf - u z ) - a z (vg - fif)g 

Another formulation of nonlinear superposition principle brings to a discrete 4-point 

(i k) 

equation on the square grid for some new variable z„ (the subscript corresponds to the 
shift in the Volterra lattice and is dummy, superscripts enumerate the Darboux transfor- 
mations). This equation is not too convenient for the purpose of the vector generalizations 
which we have in mind, however it is of interest by itself and we spend some space to 
describe it. The form of the equation depends on the sign of a. 
In the simplest case a = equations (J7]) imply the relation 

JL - JL. - f - f 

— Jn+l Jn—1 

Vn Vn-l 

which allows to introduce the variable z n accordingly to the equations 

f n = H(z n - Zn-\), 1/V n = Z n+1 - Z n -\. 

This change turns the relations ([7]) into a single equation 

{z n +l — Z n ){z n j r l — Z n ) = \l 2 
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which define Darboux transformation in terms of the variable z. Now, consider another 
Darboux transformation corresponding to the value A = v: 

{z n +l — z n){ z n+l — Z n ) = V 2 . 

The easy calculation proves that the double Darboux transformations coincide: z n = z n 
and moreover, the common value is given by the superposition formula: 

{z n z n ){z n Z r i) — /i V 

In other words, the Darboux transformations and superposition formula form the triple 
which is 3D-consistent, or consistent around a cube [57]. The iterations of Darboux trans- 
formation bring to the discrete equation on the square grid (with fixed subscript) 

(z o-+i, fe+ D _ Z m ){z u,k + D _ z u + i, k)) = ( ^))- 2 _ (i/W) - 2 . 

This is a very well-known 3D-consistent equation which defines as well the nonlinear su- 
perposition principle of the classical Darboux transformation for Schrodinger operator. 
This coincidence is not too surprising since it is known for long that Volterra type lat- 
tices are symmetries of the dressing chains which define Backlund transformations for 
KdV type equations (this relation was discussed, from the different points of view, e.g. in 

pa [annum). 

Analogously, in the case a = — c 2 the variable z n is introduced accordingly to the 
formulae 

_ fl(z n + Z n -\) _ Z n+ i + Z n -x 

Jn — / ~ \ i v n — C 

C\Z n Z n —\) Z n -\-± Z n —\ 

After this the relations ([7J turn into equation 

a (z n +l — Z n )(z n+ i — Zn) = — fl 2 (z n +l + Z n )(z n +1 + Z n ), 

and (fTTj) is replaced by equation 

r{^){z^z^ + ^•+ 1 '*)^+ 1 ' fc + 1 )) = r(u^){z^ k h^ + Z V> k+X > zU +1 > k+ V) (12) 

where r(A) = (A 2 — a)/(A 2 + a), which is equivalent to nonlinear superposition principle 
for sinh- Gordon equation. 

Finally, if a = c 2 then the change 

, M(l + Zn— l z n ) 1 + Zn+lZn— 1 

In — 7 _ x > V n — C 

C\Zn— 1 Z n ) Z n -\-l Z n — 1 

is used which brings equations ([7]) to the form 

a{z n+ i - z n ){z n +i ~ z n ) = fJ, 2 (l + z n+ iz n )(l + z n+ iz n ) 
and leads to equation 

(r(M«) + r{u^)) - zg+W) (z^ - zg+W) 

= (r(MW) - r(*f>)) (1 + 4 J ' fc+1) 4 J+1 ' fc) ) (1 + 4' fc) 4 J+1 ' fc+1) ) 

equivalent to nonlinear superposition principle for sin-Gordon equation. Equation (\12\) 
turns into (fT3l) under the complex change z — ► (i — ^)/(i + z), so that these equations are 
two different real forms of one and the same equation. 
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Concluding this Section, we notice that an analogous construction scheme exists also 
for solutions of the Volterra lattice 

U n ,x = U n (u n+ l - U n -l). 

The corresponding formulae are even much simpler, for example the equations 

U n = (v n ~ fl>)(v n +i +fi), U n = (v n+ i - fl)(v n + (j) 

replaces ((ZJ) while the role of the lattice © is played by the lattice flU) at a = —/J- 2 . 
Therefore, the lattice Q is actually the second modification of Volterra lattice. This 
sequence is analogous to the sequence of equations KdV — > mKdV — > exp-CD which can 
be obtained by continuous limit from the lattices under consideration. Unfortunately, 
although the Volterra lattice admits some multi-component generalizations [28, 29J, the 
vector ones are absent, this is why we have started from the more complicated object. 



3 First vector generalization 



Sometimes the zero curvature representation for a vector generalization can be obtained 
just by passing to the block matrices. Unfortunately, this is not the case for the matrices 
([5]), (jSJ). It turns out, however, that such block generalization is easy if one consider the 
linear equations for the length 3 vector with the components consisting of the products of 
the components of ^f. Additionally, it is convenient to apply a gauge transformation in 
order to make the determinants of the matrices L, M constant and the matrix A traceless. 
In this way we come to the following matrices which define, as can be easily verified, the 
zero curvature representation for the lattice at a = and its Backlund transformation: 



L, 





v 1 






1\ 




A 


-1 








2A 


X 2 


V n 


vl) 



I A 2 



A 



An = A 



V 

/ -A 2v n -i 

-V n 



A*/r7 
1 



2A 
-1 — — 

2Xf n 



1 

Xfn 



\ 



o \ 

V n ~l 



\ -2v n A J 

It is not a problem to find the matrices for the general case a ^ 0, but they are more 
cumbersome. Fortunately, we will not need them, since one of the vector lattices exists 
only in the case a = anyway, and for the second one this assumption does not lead to 
the loss of generality (see Section . 

The block matrices for the vector lattices are derived from here under the "proper" 
interpretation of vector-valued quantity. To make notation more clear we toggle 

to the upper case for the vectors. We assume that the vector space is equipped with a 
symmetric scalar product (U,V) = (V,U). The identity operator is denoted I and the 
linear form U T , inverse vector V~ 1 and operator UV T are defined as follows: 

V T (U) = (V,U), V~ 1 = 1 y— UV T (W) = U{V,W). 
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In the case of finite-dimensional Euclidean space one can think of V as of the column 
vector and of V T as of the row vector. 

The first vector analog of the lattice (j3|) exists only at a = 0. It is of the form [S] 



V, 



TLX 



2(K, V n+ i - V n - X )V n - (V n , V n )(V n+ i - V n -i). 



This lattice appears as the compatibility condition for the linear systems 



D, 



Vn-1 
^n 



1 

-I xv- 1 

-2A(F- 1 ) T \ 2 /{V n ,V r 




-A 2V n T _ 1 

-V n V n -i 

-2Vj A 



Vn-1 
*n 



(14) 



(15) 



(16) 



Notice that the systems (|15p . (|16|) possess the first integral in common 
J = (* n> * n ) - ^V„-i, (T - 1)(J) = D«(J) = 0. 



(17) 



The Darboux transformation is defined by a particular solution at the zero level of this 
first integral. 

Statement 1. Let F n = & n /(f) n where ip = <ft, \£ = <!> is a particular solution of the linear 
systems \15\). \16\) at A = \x and at J = 0. Then the transformation 



( 




A 2 



V 2 (F n ,F n ) 

V 1 



A 2 



V 



1 



H { n ' 

9X 2 
,,2 n n 



1 



\ 



A 2 

o (F n , -^n) 
^ 2 




(18) 



maps £/ie general solution of these systems into the solution of the systems of the same 
form, with the original and transformed potentials related by equations 



fjV n 1 — F n + F n+1 , fxV n 1 — F n+ i + F n 
The expanded form of relations (|19|) is 

F n +1 + (F n+ i, -F n+ i)F n 



-1 



(19) 



V n = f-L 
V n = fJL 



1 + 2(F n ,F n+ i) + (F n , F n ) (F n+ i, F n+ i) 
F n + {F n , F n )F n+ \ 



(20) 



1 + 2(F n , F n +-\) + (F n ,F n )(F n 4-i,F n +i) 
Each of these transformations is the substitution to the lattice (|14h from the lattice 

F n , x = /i 2 (Fn-i + F,- 1 )- 1 - fi 2 (F n+1 + F" 1 )- 1 . 
It is easy to see that these formulae turn into (J7|) and ([9]) in the scalar case at a = 0. 
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Figure 2. A soliton of the lattice JUJ); C = (1,0), If = (0,1), 7 = 2, 

ci = 1, c 3 = -1. 



The derivation of the nonlinear superposition principle is not more difficult than in the 
scalar case. One comes by multiplying the matrices M of the form (|18p to the following 
Yang-Baxter mapping (cf eq. (fTUj) at a = 0) 

mrr , (/z 2 - ^ 2 )(G, G)F + F) - 2 j »(F, g) + G))G 
R V?,G,hu) = (G,G)(vF-»G,vF-pG) ' (21) 

It is possible to obtain the analog of equation (jllj) . too. Let us introduce the new 
vector variable Z n accordingly to the formulae 

F n = ^(Z n — Z n -{), V' 1 = Z n+ i — Z n -\. 

Equations (|19|) become equivalent to the single equation 

Z n +i — z n = n 2 (Z n+ i — Z n ) 1 

under this change. Next, consider Darboux transformation corresponding to the spectral 
value A = v. 

Z n +i — Z n = v 2 {Z n+ \ — Z n ) 1 . 

The direct calculation shows that the repeated Darboux transformations coincide: Z n = Z n 
and moreover, the result is given by the equation 

Z n ~ Z n = (/i 2 — v 2 )(Z n — Z n ) 1 . 
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Figure 3. A breather of the lattice (HH); C = (1,0), K = (0, 1), 7W = 
7^=0.5 + *, «4 1 )=ci a) =0.5 J c^cf^l. 

Iterations of the Darboux transformation are governed by the 3D-consistent discrete equa- 
tion on the square grid (the subscript n is dummy): 

z u+i,k+i) _ z m = ( (/U (i))-2 _ {y{k)y^ z ij,k+i) _ Z ii+my\ 

This equation with important applications in the discrete geometry was introduced in [12] 
(a special reduction was considered in [Taj ) , see also [T3J [J3] . 

Let us make use of Darboux transformation for construction of the soliton solution. 
The solution of the linear equations (|15p . (|16p with constant coefficients V n = C = const, 
(C, C) = 1, at A = /i reads 

^ = Cl7 " e (7-7- 1 )- + C27 -« e (7- 1 -7)- + 2C3 , 
$ n = (-l) n K + 11 e (7-7- 1 )- + ° 27 " e (7- 1 "7)- + C3 ^j C (22) 

\l + 7 1 + 7 1 / 

/U 2 = 7 + 2 + 7 -i ) (C,K)=0, 7 (K,K) = (l- 7 ) 2 (cic 2 -c|) 

(the latter relation is equivalent to the constraint J = 0; and we do not consider the 
cases of multiple eigenvalues fi = 0, /x = ±2). The equations ([19]) bring, after elementary 
transformations, to the one-soliton solution of the lattice (fig. [2]) 

n — J.2 

Clearly, this solution always lies in the plane of the vectors C, K, that is it is actually 
2-component, independently on the dimension of the vector space under consideration. 
The C-component is a soliton on the unit background. Its shape is slightly different for 
positive and negative values of C3. The i^-component has localized oscillations on the zero 
background. They originate from the powers of —1 in the solution (|22p rather than a pair 
of complex conjugated points of the discrete spectrum, that is this solution is not a genuine 
breather. However, the additional dimension makes the breathers possible as well, in spite 
of the absence of the parameter a (fig. [3]). The iV-soliton solution is constructed starting 
from a set of solutions of the form (|22p and evolves in the space spanned over the vectors 



4 Second vector generalization 



The lattice 

U n , x = ((U n , U n ) + a)(U n+1 - Un-i) (23) 

looks more natural and simple generalization of the lattice In contrast with (|14p it is 
integrable at arbitrary value of parameter a, though it turns out to be not so important as 
in the scalar case. Indeed, it can be easily eliminated or, more precisely, "confined inside 
the lattice" at the expense of increasing by 1 the dimension of the vector space under 
consideration. This is done by means of the orthogonal complement: let U n be a solution 
of the lattice (1231) then the vector 



V n = U n + E, E = const, (U n ,E) = 0, (E,E)=a 
(if a < then a pseudoeuclidean scalar product is used) satisfies the lattice 

Vn.x = {V n ,V n )(V n +i — V n -i). 



(24) 



(25) 



This transformation does not lead to any problem when constructing solutions since the 
reduction (|24l) is consistent with higher symmetries and Backlund transformation. On the 
other hand, all formulae simplify essentially (cf e.g. equations ([29]) and ([33]) below). The 
matrices of the zero curvature representation become simpler as well. 
The lattice (I25D is the compatibility condition of the linear systems 




/0 

I-2V- 1 VJ 



i 



\ 



x 2 /{v n ,v n )J 




(26) 



Dr 



Vn-1 



2AK 



-\V n 2V n V n T _ 1 



n-l 



2K-iK T 






\ 










XVn- 


-1 




X 2 






) 



(27) 



V -2AF n T 
These systems possess the first integral (|17p in common, like in the previous case. 



Statement 2. Let F n = & n /4> n where ifj = (f),ty = <&isa particular solution of the linear 
systems HZo)) . jIT7[) at A = \x and such that J = $ n ) — 4> n 4>n~i = 0. Then the transform 



( 




A 2 



V 2 (F n ,F n ) 



2A 
A* 



\ 



[i 2 



2F~ l F 



\ 



1 



2 A F T 



-F n 

A 2 

9 {Fni F n ) 

LI 2 




(28) 



maps the general solution of these systems into solution of the systems of the same kind, 
with the original and transformed potential related by the equations 



Vn 



F n +i — (F n+ i,F n+ i)F n 
1 — (F n ,F n )(F n+ i,F n+ i)' 



V n = fi 



Fn {FniF n )F n +l 

1 — (F n ,F n )(F n+ i,F n +i) 



(29) 
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In comparison with the previous Section, the equations (|29p are slightly shorter than 
20]) . but the analog of the lattice © is more cumbersome: 

p n 2 (F n , F n )((F n — -P^i) -1 - (F„ -F-^)- 1 ) 



The superposition of Darboux transformations is defined by the Yang-Baxter map 

F ti,k,«) = R ^),F^;^\^), = R(F^,F^;^ k \^), 

R i F > G >P> v ) = (G,G)(vF-»G,vF-pG) ' (30) 



Notice that the formulae (|21jl and (I30p coincide up to the permutation of \x and v in the 
numerator. Despite of such similarity, an analog of equation (jlip is probably lacked in this 
case. 

Statement 3. The Darboux transformation is consistent with the reduction \21$ . 

Proof. Let us apply the change V n = U n + E, F n = H n + h n E, (H n , E) = to the equations 
(|29p . with the scalar factor h n unknown for the moment: 

TT , _ H n+ i + h n+ \E — ((H n+ i, H n+ \) + ah n+1 )(H n + h n E) 
U n + E = fi- 



U n + E = fi 



1 - ((H n , H n ) + ah 2 )((H n+1 ,H n+1 ) + ah? n+l ) 
H n + h n E - ((Hn, H n ) + ah 2 n )(H n+1 + h n+1 E) 
1 - ((H n ,H n ) + ah n )((H n+1 ,H n+1 ) + ah 2 n+1 ) ' 



(31) 



Collecting the coefficients of F yields the coupled algebraic equations for h n and h n +\. 
It is not obvious beforehand that their solution is compatible with the shift in n. If this 
would be not the case then the change F n = H n + h n E would be incorrect. However, the 
direct computation proves that h n is defined by one and the same formula for all n as a 
solution of quadratic equation 

ah 2 n - fih n + (H n , H n ) + 1 = 0, (32) 

therefore the ^'-component is detached in the transformation (I29p . □ 

Statement [3] makes abundant the separate study of the case a / 0. Nevertheless, 
all formulae can be, in principle, rewritten for this well, moreover, their rational 

structure can be preserved by use of the stereographic projection for the quadric ([320 : 

(v 2 - a)F n v(l + (F n ,F n )) _ a 

Mn ~ u 2 + a(F n ,F n y z/ 2 + a(F n , F n ) ^~ V+ u 

(more rigorously, some other letter should be used here instead of F, but we hope it will 
not lead to misunderstanding). For instance, the substitution into (|3ip brings, under this 
parametrization, to the Backlund transformation for the lattice with parameter ([23]): 

= (v 2 + a(F n , F n ))F n+1 - (a + v 2 (F n+1 , F n+1 ))F n 
u(l-(F n ,F n )(F n+1 ,F n+1 )) ' 
~ = (v 2 + a(F n+1 ,F n+1 ))F n - (a + v 2 (F ni F n ))F n+1 
u(l-(F n ,F n )(F n+1 ,F n+1 )) 
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The Yang-Baxter map fl3Qf) can be rewritten in more general form in a similar way. The 
transformation ()33[) turns into f|29|) at a = 0, while in the scalar case we come back to the 
transformation (|7|), under identifying U,F,v with v,f,fj, respectively. 
It is not difficult to compute, by use of ([29]) . the one-soliton solution 

where the solution (f> n , <3? n of equations ([2^1 , (|27p with constant coefficients is given by 
almost the same formulae (122 H as before, with the only difference that the factor (— l) n in 
front of K disappears. This distinction leads to the absence of blinking oscillations in the 
one-soliton solution which is more natural from the point of view of the continuous limit. 
In all other respects the construction of multisoliton and breather solutions is analogous 
to the previous case. 



5 Higher symmetries and associated systems 

Both vector lattices (|14jl and (|25p belong to an infinite hierarchy of commuting flows. 
We restrict ourselves by consideration of the simplest higher symmetries which are of the 
second order with respect to the shift in n. In the scalar case one has, setting a = for 
simplicity, the pair of consistent lattices 

Vn,x = Vn(v n +l-V n -i), V n ,t = Vn(l% +1 (v n +2 + «n) ~ ^-1 ( v n + v n-2))- (34) 
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Obviously, the first of these equations can be solved with respect to v n+ \ or i? n _i, and this 
allows to express recursively all Vj through the pair of variables u = i> n +i> v = v n . After 
this, the symmetry takes the form of Kaup-Newell evolution system 

u t = u xx + (2u 2 v) x , v t = -v xx + (2uv 2 ) x 

and the shift in n defines an explicit auto-substitution for this system (the simplest type 
of Backlund transforms). 

The lattice (|14p and its symmetry can be compactly written in the form preserving the 
structure (j3l1) 

V n , x = P Vn (V n+1 - V n -t), V n , t = P Vn (P Vn+1 (V n +2 + V n ) - iV n _! (V n + V n - 2 )) (35) 

by use of the operator Py(U) = 2(V, U)V — (V, V)U. It is easy to check that the identity 
(Py)~ l = Py-i is valid for this operator. Making use of it one can solve, like before, 
the first equation with respect to V n+ \ or V n -\ and to express all Vj through the pair 
of variables U = V n +\, V = V n . This brings the symmetry to the form of the vector 
generalization of Kaup-Newell system 

U t = U xx + (4(U,V)U-2(U,U)V) x , V t = -V xx + (A(U,V)V -2(V,V)U) X . (36) 

Analogously, the commuting flows for the second vector lattice are 

V n>x = (V n ,V n )(V n+1 -V n ^), 

V ntt = (V n , V n ) {(V n+ i, V n+1 )(V n+2 - V n ) + (K-i, V n -i)(V n - V n - 2 ) (37) 
+ 2((V n+ i,V n ) + (V n , - K-l)) 

and the associated evolution system reads 

U t = U xx + A(U,V)U x + 2(U,U)V x , V t = -V xx +4(U,V)V X + 2(V,V)U X (38) 

which is another vector analog of Kaup-Newell system. 

Another interesting type of associated systems is obtained for the scalar quantities 

Vn = (V n , V n ), q n = 2{V n , V n -i) 

which satisfy, in virtue of any of the pair (|35H or ()37l) one and the same two-dimensional 
modified Volterra lattice 

Pn,t + 2p n (p n+ i - p n -l) = p n ((q n +l + q n )x + q n +l ~ qI), Pn,x = Pn(q n +1 ~ ?n)- (39) 
It can be written in the form 

Pn,t + 2p n (p n+ l - Pn-l) = (r n Pn)x, (Pn+lPn)x = Pn+lPn(r n +l ~ r n ), 

as well, where r n = q n +i + q n = 2(V n , V n+ \ + V n -i)- These lattices are closely related to 
Mikhailov lattices introduced in |18j . 

The lattices (|35l) . (f37|) can be effectively used for the construction of particular solutions 
of the systems ([36j) . (f38|) and the lattice (f39j) . Along with the construction method of 
the soliton-type solutions described above, one can use to this end the periodic closure 
V n +N = CV n with orthogonal operator C which leads to the finite-dimensional dynamical 
systems. 
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6 Further vector analogs 



Remind that the classification problem of scalar integrable lattices of Volterra type was 
solved by Yamilov j2U] within the symmetry approach. Recently one of the authors has 
obtained an analogous classification of the vector Volterra lattices on the sphere, that is 
under the constraint (V n ,V n ) = 1 [25]. This constraint essentially simplifies the problem 
which is very complicated and remains open for the case of free space. Other simplifying 
assumptions can be used of course, for example the polynomiality of the lattice. It should 
be noted that in the continuous case very many polynomial equations are known; we 
mention only the papers [H0E] which contain the examples and some classification results 
for the vector systems of derivative nonlinear Schrddinger type, equations (|36|) . (|38l) being 
just two instances of such systems. In the discrete case, however, the polynomiality is not 
too natural assumption, as one can see already from the Yamilov list of scalar lattices. 
In the vector setting we have not succeeded in finding another polynomial Volterra type 
lattices possessing higher symmetries aside from (|14p . (f23"j) . 

Our search of integrable lattices was based on the straightforward method of undeter- 
mined coefficients. In the simplest case the lattice and its symmetry are of the form 

V n , x = aVv n+l + a^V n + a^V n - U V n , t = 6 {2 V„ +2 + • • • + b^V n . 2 

where the scalar coefficients a® are linear with respect to the scalar products of V n +i, V n , 
V n -\ and &W are quadratic with respect to the scalar products of V n+ 2, ■ ■ ■ , V n -2- It is easy 
to find that the homogeneous lattice contains 18 parameters and its symmetry contains 
600 ones. Calculating of the cross derivatives yields a system of bilinear equations for the 
coefficients. Although this system is very bulky, its solving is, in principle, not difficult 
since the equations are very overdetermined and sparse (in particular, a large part of 
equations is monomial). The answer is the consistent pairs of the lattices (|35|) and (|37|) 
(there are also few solutions with a' 1 ' = a^ 1 ) = 0, but all such lattices can be reduced to 
the scalar ones and therefore they are not of interest for us). 

It is clear that the scope of this method in this problem is very restricted. If one takes 
the coefficients a" quadratic with respect to the scalar products and foW of the fourth 
degree then the number of unknown parameters in the lattice and its symmetry becomes 
63 and 15300 respectively, and even the calculation of the commutator becomes not so 
trivial task. This case is still manageable, but with the empty answer. 

We also have partially analyzed the case when the lattice is of the second order with 
respect to the shift in n and its symmetry is of the fourth order, that is 

Vn,x = a (2 V n+2 + • • • + a(" 2 V n -2, Vn.t = fr (4 V n+4 + • • • + fr ( " 4) K_ 4 . 

One may hope that some vector analogs of Narita-Bogoyavlensky lattice [301 [3T| [29] appear 
here, more precisely, analogs of some its modification with odd degree of nonlinearity, for 
example 

V n ,x = V n (v n+ 2V n+ l - V n -lV n ^ 2 ) Or V UjX = V n+ lV^V n -.l{v n+2 V n+ l - V n -lV n -2) ■ 

Notice that classification of such lattices is not known even in the scalar case. Unfortu- 
nately, the analogs of Narita-Bogoyavlensky lattice have not been discovered, however we 
have found two more lattices relative to Volterra lattice: 

v n , x = (v n ,v n )((v n+1 ,v n+1 )(v n+ 2 + v n ) - (v n -i,v n -i)(y n + v n - 2 )), (40) 

V n , x = {V n+1 , V n )(V n , " Vn-2). (41) 
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Each of these lattices possesses 4-th order symmetry which we do not bring because of their 
length. The study of these examples falls beyond the scope of our article. We only notice 
that the lattice (|4(J[) generalizes the second flow of the modified Volterra lattice (|34p . so 
that this flow admits at least three vector analogs. The question on the number of vector 
analogs for the higher flows of the hierarchy remains open. The lattice (|41|) in the scalar 
case is a modification of the Volterra lattice on the "stretched" grid: 

V n ,x = V n+ iV n V n -i{V n -2 ~ Vn-2) > V>n,x = U n (U n -2 - U n _ 2 ) , 

but in the vector case this substitution makes no sense and the lattice (I4ip seems to be 
an independent object. The zero curvature representations and Backlund transformations 
for the lattices ([4*01) . (l4*Tj) are not known for now. 

Summing up, we may say that the classification of the polynomial lattices of Volterra 
and Narita-Bogoyavlensky types is a very difficult open problem, probably with very scarce 
answers. The alternative approaches to the method of undetermined coefficients are the 
analysis of the necessary integrability conditions in the form of canonical conservation 
laws [201 IHl E2] and the perturbative approach [32J, however the contemporary state of 
the theory does not allow to effectively apply them, even in the scalar case. 
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